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Abstract With the discovery of Dark Energy, Ade, there is now a universal 
length scale, £de = c/(AdeG) 1 ^ 2 : associated with the universe that allows for 
an extension of the geodesic equations of motion. In this paper, we will study 
a specific class of such extensions, and show that contrary to expectations, 
they are not automatically ruled out by either theoretical considerations or 
ON ' experimental constraints. In particular, we show that while these extensions 

affect the motion of massive particles, the motion of massless particles are not 
changed; such phenomena as gravitational lensing remain unchanged. We also 
show that these extensions do not violate the equivalence principal, and that 
because £de — 14010§2o Mpc, a specific choice of this extension can be made 
so that effects of this extension are not be measurable either from terrestrial 
experiments, or through observations of the motion of solar system bodies. A 
lower bound for the only parameter used in this extension is set. 
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OQ ' 1 Introduction 

in 

(**»» ■ The recent discovery of Dark Energy (see pQ, [2] and references therein) has 

broadened our knowledge of the universe, and has demonstrated once again 
that it can hold surprises for us. This discovery has, most assuredly, also 
brought into sharp relief the degree of our understanding of the universe. In 
this paper, we will study one specific implication of this discovery: With the 
discovery of Dark Energy, Aue , there is now a universal length scale, £de = 
cZ(AdeG) 1 / 2 , Q|,[j associated with the universe that allows for extensions of the 
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1 Our £oe differs from the length scale, r^, defined in [3] by a factor of %/3- 

2 ft is also possible to construct the length scale (Hc/A^e) 1 ' 4 ' ~ 85 fi m. Experiments 
have shown that this scale does not affect the Newtonian potential [!]• 



geodesic equations of motion (GEOM). We find that contrary to expectations, 
these extensions are not automatically ruled out by theoretical considerations, 
nor are they ruled out by experimental constraints either through terrestrial 
experiments or through solar system tests of general relativity. Indeed, we show 
in this paper that one specific extension of the GEOM is a viable alternative 
to the GEOM, and we obtain a lower bound for the only free parameter used 
in its construction, a power-law exponent, cxa- 

There are good theoretical and physical reasons for studying the range of 
extensions of the geodesic equations of motion that are allowed. Arguments 
for the use of the geodesic equations of motion to describe the motion of 
massive test particles in curved spacetime are based on various statements of 
the equivalence principle, and the principle of general covariance (see chapter 4 
of [5]), along with arguments in favor of simplicity and aesthetics. Importantly, 
these arguments are made in addition to those made in favor of Einstein's field 
equation. Namely, these is no unique way of deriving the geodesic equations of 
motion from the field equations. Indeed, in 1938 Einstein, Infield and Hoffman 
attempted to show that as a consequence of the field equations, massive test 
particles will travel along geodesies in the spacetime [B] . These attempts have 
continued to the present day [7], [8]. 

Extensions and modifications of the GEOM have been made before, of 
course. On the level of Newtonian dynamics, Modified Newtonian Dynamics 
(MOND) [5] has been proposed as an alternate explanation of the galactic 
rotation curves. On the relativistic level, there has been recent efforts [TU] to 
develop a general framework to study modifications to the GEOM in the weak 
field, linearized gravity limits. The major impetus for this work has been to 
describe a series of dynamical anomalies — the Pioneer anomaly (see [TT] and 
|12j). the flyby anomaly [13], and the lengthening of the Astronomical Unit 
[2] — that have been observed at the Solar system scale. 

The focus of this paper is to establish the underlying theoretical framework 
that can be used to describe structures and dynamics at the galactic length 
scale and above. In a future paper [15], this framework will be applied to an 
analysis of the galactic rotation curves, and the impact that this extension 
has on phenomena at cosmological length scales will be studied. As such, we 
focus here on the Dark Energy energy length scale, and on how the existence 
of this scale allows for extensions of the GEOM. Indeed, we find that with 
this length scale, £de, extensions of the GEOM are not difficult to construct. 
The quotient c 2 R/ AdeG is dimcnsionless, and functionals of this quotient can 
easily be used to extend the GEOM. What is more relevant is whether or not 
the resultant equations of motion will be a physically viable alternative to 
the GEOM. As such we will be guided in our extension of the GEOM by the 
four conditions listed below. They are deliberately chosen to be conservative 
in scope, and thus stringent in their application. Somewhat surprisingly, we 
will show that there is at least one extension of the GEOM that satisfies all 
four. 

First, we require that the extension preserve the equivalence principle, 
which is one of the underlying principles on which general relativity is founded. 



In the following sections of the paper, we will explicitly see that this preser- 
vation is assured by the fact that Ide is the same for all test particles. This 
universal nature of Ide is crucial. While other length scales — say, the pro- 
ton mass — could be used for the extension, the resultant equations of motion 
would depend on this mass. They could not then be applied to the motion of 
protons without explicitly violating the uniqueness of free fall condition. 

Second, we require that the extension not change the equations of motion 
for massless test particles; such particles must still follow the GEOM. All astro- 
nomical observations — of which gravitational lensing is playing an increasingly 
important role — are based on the motion of photons of various wavelengths. 
Modifications to the equations of motion for photons will require a reintcr- 
pretation of these observations, a daunting step not to be taken without good 
reason. We will show that by considering a class of extensions that is based 
on conformally scaling the rest mass of the test particle, we arrive at extended 
GEOMs that, on the one hand, will not change the motion of massless test 
particles, but will, on the other, change the motion of massive ones. Being con- 
formal, the motion of photons will not be affected by this class of extensions, 
and they will still travel along null geodesies; phenomena such as gravitational 
lensing will remain unchanged. While in form this class of extensions resem- 
bles a scalar field theory that is non-minimally and nonlinearly coupled to the 
scalar curvature, R, such theories are constructed at the quantum level. Our 
extension of the GEOM is done at the classical, h = 0, level, with the scale of 
the coupling set by Ide- 

The third condition involves the attempts [5], [7J, [5] at proving that the 
GEOM are the unique consequence of the Einstein field equations (see also 
page 72 of [5]). These proofs would seem to rule out any physically relevant 
extension of the GEOM, and by necessity, our extension of the GEOM cannot 
be precluded by such proofs. That the extension is possible is because these 
proofs focus on the motion of test particles in regions where the Einstein tensor, 
G^y, vanishes. We will see that in these regions the extended GEOM reduce 
to the GEOM, and thus do not violate these proofs. Indeed, we will explicitly 
construct the energy- momentum tensor for an inviscid fluid of massive particles 
propagating under the extended GEOM. 

The fourth condition is the most stringent of the four. With the excep- 
tion of the as- yet unexplained anomalies described above, we require that the 
extension of the GEOM not produce effects that are measurable either in ter- 
restrial experiments, or through the motion of bodies in the solar system that 
have traditionally been used to test general relativity. While stringent, we will 
nevertheless show explicitly that a choice of extension can be made which satis- 
fies it. Physically, this choice is possible because at (7.2li ' g4 ) x 10 -30 g/cm 3 
[16], Ade is orders of magnitude smaller than the density of matter, pihmt, 
cither currently achievable in terrestrial experiments (where densities exceed 
10~ 18 g/cm 3 ), or present in the solar system (where the density of matter 
in Mercury's orbit exceeds 10 -23 g/cm 3 ). Correspondingly, at MOlOg™ Mpc 
Ide is more than three times larger than the observed size of the universe, 
and is orders of magnitude larger than the solar system. Nevertheless, we find 



that even though the disparity between the magnitude of Aoe and pii m it — or, 
equivalently, between ioE and the size of the solar system — is large, a non- 
linear function of c 2 R/ AeeG is needed in constructing the extension for its 
effects not to have already been seen in terrestrial experiments. The simplest of 
these extensions has only one free parameter, a a, a power-law exponent that 
determines the behavior of the function at densities both much larger than 
Ade, and much smaller than it. Lower bounds for a a are set by requiring 
that the extension does not produce observable effects in current terrestrial 
experiments. 

While it may be possible to apply the analysis in this paper to the expla- 
nation of Solar system anomalies such as the Pioneer anomaly, the focus of 
this research is on phenomenon at the galactic scale or longer. It is for this 
reason that we require our extension to be constrained only by experiments 
and observations that are currently well-understood, and for which the under- 
lying physics is well-known. We leave to future work the question of whether 
or not our analysis can be applied to explaining the Pioneer and other Solar- 
system-scale anomalies. 



2 Extending the Geodesic Lagrangian 

We begin our extension of the GEOM with Einstein's field equation in the 
presence of a cosmological constant 

p _I p i AdeG _ ^G ^ 

where T^ is the energy-momentum tensor for matter, R^ v is the Ricci tensor, 
Greek indices run from to 3, and the signature of g^ is (1, — 1, — 1, —1). 
While there is currently no consensus as to the nature of Dark Energy (pro- 
posals have been made that identify it with the cosmological constant Aee 
[3] , with quintessence [T7] , [TB] , [T2] , or even as a consequence of loop quantum 
cosmology [20]), modifications to Einstein's equations to include the cosmo- 
logical constant are well known and are minimal. We will thus identify Dark 
Energy with the cosmological constant in this paper, and require only that 
Ae>e changes so slowly that it can be considered a constant in our analysis. 

Requiring that Eq. (fT]) still holds under the extension of the GEOM is a 
choice, one which, we will see below, is the simplest. Although it may seem 
surprising that we can still make this choice even though we will be changing 
the GEOM, extensions of the GEOM need not change the relation between 
R^ v and !),„ given in Eq. (TT|). They can rather change the precise form that 
T^ v takes for matter. To see this, consider the following. 

The total action, S, for a system consisting of gravity, radiation, and mat- 
ter can be written as a sum of three parts: S — Sgrav + Radiation + ^matter- 
Here, Sgrav is the action for gravity, SVadiation is the action for radiation, and 
Smatter is the action for matter. We will show below that the class of ex- 
tended GEOM we consider here will not change the equations of motion for 



radiation so that ^radiation w ill not be changed. The extension will certainly 
change Scatter, however, and in the next section we will explicitly construct 
the energy-momentum tensor for an inviscid fluid whose constituents follow 
an extended GEOM. For the present argument, we only need to note that 
whatever the form taken for ^matter, we are still free to choose S grav to be the 
Hilbert action; Eq. ([T]) then follows after taking the functional derivative of S 
with respect to the metric. 

This ability to change Scatter while leaving S gra _ v , and thus Eq. (p}, un- 
changed was explicitly exploited in the construction of minimally coupled 
scalar fields, 4>r- There, the mass term of the scalar field, m 2 0|j, is replaced by 
Rm 2 (f) 2 R in the action for matter, and yet S^rav is still taken to be the Hilbert 
action (see Sec. II. D). Einstein's field equations, Eq. ([T]), still hold; the only 
change is the form that T^ takes. 

Both the geodesic Lagrangian 

dx» dx v x 1/2 



and the GEOM 



£o = mc(g» u — — ) , (2) 



v v V v v» = — = 0, (3) 

(where v^ = x^ is the four- velocity of the test particle), have natural geometric 
meaning. The first is a proper time interval, while the second is the equation 
for parallel transport, which determines the shortest time-like path connecting 
two points. But, aside from their inherent geometric meaning, there is also a 
good physical reason to take Eq. ([3]) as the equations of motion for a test 
particle. In the absence of Dark Energy, Eq. © is the most general form that 
a second-order evolution equation for a test particle can take which still obeys 
the equivalence principle. 

Any extension of Co would require a dimensionless, scalar function of some 
fundamental property of the spacetime folded in with some physical property 
of the universe. In our homogeneous and isotropic universe, there are few op- 
portunities to do this. A fundamental vector certainly does not exist in the 
spacetime, and while there is a scalar (the scalar curvature, R) and three ten- 
sors (g^v, the Riemann tensor, R^ v ^ , and the Ricci tensor, R^), R^,^ has 
units of inverse length squared. While it is possible to construct a dimension- 
less scalar m 2 G 2 R/c 4 for the test particle, augmenting £ using a function of 
this scalar would introduce additional forces that will depend on the mass of 
the test particle, and thus violate the uniqueness of free fall principle. It is 
also possible to construct the scalar g^v^v 11 /c 2 , but because of the mass-shell 
condition, v^v^ = c 2 , any such extension of £ will not change the GEOM. 
Scalars may also be constructed from R^ and powers of R„ v J^ by contract- 
ing them with the appropriate number of u M /c's, but these scalars will once 
again have dimension of inverse length raised to some power, and, as with the 
Ricci scalar, once again a rest mass m is needed to construct the dimensionless 
quantity. 



The situation changes dramatically in the presence of Dark Energy. With 
a universal length scale, £de, it is now possible to construct from the Riemann 
tensor and its contractions dimensionless scalars of the form, 

c 2 R R^v v gW Qj/37 v»v»v y v s S0 

A DE G' A DE G ' {A DB Gy K,M *^ Hv ' (4 B C) aW """ W 

Although extensions to Co can be constructed with any of these terms, we 
look at extensions with the form: 



-^Ext = 



l + V(c 2 R/A DE G) J [ 9liV ^-— ) = V\[c 2 R/A DE G]£o, (5) 



with the implicit condition that v 2 — c 2 for massive test particles. We make 
this choice for the following reasons. 

First, having 9\ be a function only of c 2 R/A EE G is the simplest extension 
that can be chosen; other choices will induce velocity-dependent effects in the 
extended GEOM. Second, we will find below that any extension of Co of the 
form Eq. ([3]) will not change the equations of motion for massless test particles; 
they can still be reduced to the GEOM. Extensions of the form Eq. (0 will 
only affect the motion of massive test particles. Third, we require that the 
extension of the GEOM not produce effects that should have already been seen 
in terrestrial experiments; these experiments are done in the nonrelativistic 
and weak gravity limits. Constraints in the choice of *H are thus found in these 
limits, where it is clear from Eq. Q that the second term reduces to the first, 
while the other terms are higher order in the curvature. We are thus left with 
c 2 R/A EE G with which to construct an extension of the GEOM. 

As £de — 14010g2Q Mpc, the question remains whether it is possible to use 
another, shorter length scale in its place to extended the GEOM; this extension 
could then be used to describe deviations from geodesic motion on shorter 
length scales. One such application would be in explaining Solar system scale 
anomalies such as the Pioneer and flyby anomalies. At this scale, a natural 
length scale would be MqG/c 2 = 1 480 m, where M Q is the mass of the Sun, 
and the resultant extension of the GEOM may be applicable to the description 
of motion on the Solar system scale (see [10]). Its application to the description 
of motion at the galactic scale or longer is more problematic, however, and it 
is precisely on these length scales that we are concerned with here. On the 
galactic scale, stars can be treated as test particles, and as MqG/c 2 depends 
explicitly on the mass of the Sun, the use of this length scale in extending 
the GEOM would mean that the motion of stars in galaxies would depend on 
the mass of the Sun. This would not be physically reasonable, and would also 
violate the uniqueness of free fall condition. 

If T>(x) is the constant function, then £Ext differs from Co by an overall 
constant that can be absorbed through a reparametrization of time. Only non- 
constant 2)(x) are relevant. It is how fast 2)(x) changes that will determine its 
effect on the equations of motion, and not its overall scale. Indeed, in extending 
Cq we have in effect performed a conformal scaling of Co by replacing the 



constant rest mass m of the test particle with a curvature-dependent rest 
mass md\ [c 2 R/ 'AejeG] . All dynamical effects of this extension can therefore 
be interpreted as the rest energy gained or lost by the test particle due to the 
local curvature of the spacetime. The scale of these effects is of the order of 
mc 2 /L, where L is some relevant length scale of the dynamics. The additional 
forces from C^xt are thus potentially very large. For these effects not to have 
already been seen, T>(c 2 R/AdeG) must change very slowly at current limits 
to experimental measurements. 

As mentioned above, using Einstein's field equations, Eq. (fTJ), was a choice. 
In particular, notice that because £Ext is the result of conformally scaling the 
rest mass by 9\, we may choose to instead reduce £Ext to £q through the 
conformal transformation of the metric g^ v = IH 1 ' 2 ^^. However, doing so will 
result in a Ricci tensor R^ v that is nonlinear ly related to R^ v : 

R^ = R„ v - t^V, log^H - -g^V a V a log 9* 

+iv M io g *nv,iog$n- lg^y a iogu\v a io g u\. (6) 

o o 

where V^ is the covariant derivative for g^ v . As Eq. (TTJ) holds for R^v, it 
cannot hold for R^ v . It will instead be replaced by a nonlinear relation between 
R^ v and T M „, resulting in a higher-order theory of gravity. Thus, instead of 
choosing Eq. (1) to hold for the Ricci tensor and extend the GEOM, we could 
have chosen a higher order theory of gravity from the start while preserving 
the GEOM. We would argue, however, that making this second choice would 
change the foundations of classical general relativity — a much more drastic 
step — that would also result in a less tractable theory. 



2.1 The Extended GEOM for Massive Test Particles 



For massive test particles, the extended GEOM from £Ext is 
D 2 x» o / .... v^v" 



dt 2 



c 



2 n V 



g^ 5- )VAog<R[c 2 R/A DE G}, (7) 



where we have explicitly used v 2 — c 2 . It has a canonical momentum with a 
magnitude of 



P 2 = V^ L = m 2 c 2 



1 + D(c 2 R/A D eG) 



(8) 



and the interpretation of m9\[c 2 R/ AjjeG] as an effective rest mass can be 
readily seen. What also can be readily seen is that £ext and Eq. (JT]) have lost 
the geometrical meaning that £q and Eq. $5§ have. Namely, the worldline of a 
massive test particle is not that which minimizes the proper time between two 
points; it is instead one that is either attracted to, or repelled from (depending 
on the choice of 9\), regions where the scalar curvature is extremized. 



The dynamical implications of the new terms in Eq. ((?]), along with the 
conditions under which they are relevant, can most easily be seen by noting 
that R = 4A DE G/c 2 + 8irTG/c 4 , where T = Tjf. Then %c 2 R/A DE G] = 
*H[4 + 8ttT J Adec 2 ], where the '4' comes from the dimensionality of spacetime. 
Thus, in regions of spacetime where either T^ — or when T is a constant, 
the right hand side of Eq. Q vanishes, and our extended GEOM reduces back 
to the GEOM. 



2.2 Dynamics of Massless Test Particles 

For massless test particles, we note that ^Ext is related to Cq by a scaling of 
the rest mass. This scaling can be formally interpreted as a conformal trans- 
formation of the metric, g^ v . As null geodesies are preserved under conformal 
transformations up to an affine transformation, the extended GEOM is equiv- 
alent to the GEOM, and we still have v v V v v^ = 0. 

This result can be shown explicitly in the following analysis. Because mass- 
less particles follow null geodesies, for these particles we consider extensions 
of the form 

£l Kt =nc 2 R/A DE G]£l (9) 

where we have taken 

1 dafdz" 

c ° = 2 9 ^^r^r- (10) 

As usual, the GEOM comes from £q, while the equations of motion that come 
from £^xt are 

5CI SCl\SCZd3K r7 , 



V dt 8x^ 5xv J Sx^ dt ° M v ; 

Taking now the mass-shell condition, £^ — 0, Eq. (jlip reduces to simply 

D ( dx^ 



= ^ *R—r- ■ (12) 

dt \ dt J y ' 

After the reparametrizing dt —¥ 9idt [5], we arrive at the expected result, 
v^Vj,^ = 0. Importantly, this result means that the usual general relativistic 
effects associated with photons — the gravitational redshift and the deflection 
of light — are not effected by our extension of the GEOM. 



2.3 Impact on the Equivalence Principles 

The statements |21) of the equivalence principal we are concerned with here 
are the following: 

Uniqueness of Free Fall: It is clear from Eq. (J7| that the worldline of a freely 
falling test particle under the extended GEOM does not depend on its com- 
position or structure. 



The Weak Equivalence Principle: Our extension also satisfies the weak equiv- 
alence principle to the same level of approximation as the GEOM. The weak 
equivalence principle is based on the ability to choose a frame in a neigh- 
borhood of the worldline of the test particle where J^o « 0; the Minkowski 
metric, r]^, is then a good approximation to g^ v in this neighborhood. How- 
ever, as one deviates from this world line corrections to ?7 M „ appear, and since 
a specific coordinate system has been chosen, they appear as powers of the 
Riemann tensor (or its contractions), and its derivatives (see [21] and [22] ). 
This means that the larger the curvature, the smaller the neighborhood about 
the world line where rj^ v is a good approximation of the metric. Consequently, 
the weak equivalence principle holds up to terms first order in the curvature. 
As the additional terms in Eq. ((?]) are first order in R as well, our extension 
of the GEOM satisfies the weak equivalence principle to the same order of 
approximation as the GEOM does. 

The Strong Equivalence Principle: Because we only change the geodesic La- 
grangian, all nongravitational forces in our theory will have the same form as 
their special relativistic counterparts. Moreover, the extended GEOM reduces 
to the GEOM in the R -> limit. 



2.4 Connections with Other Theories 

As unusual as the extended GEOM, Eq. ©, may appear to be, there are 
connections between this extension and other theories. 



2-4-1 The Class of Scalar Field Theories in Curved Spacetimes 
The Klein-Gordon equation corresponding to the extended GEOM is 

c 2 R 



V 2 c 



r, 2 



1 + 2) 



AdeG 



= 0. (13) 



This is the equation of motion for a scalar field <f> that is non-minimally and 
nonlinearly coupled to R. Scalar field theories of this class have been studied 
before, the most notable of which is 

V 2 0r+(^+£R)<^ = O. (14) 

When £ = 1/6, the scalar field will be conformally invariant even though mj^O 
[23]. 

There are important differences between these theories and the theory we 
are considering here, however. Scalar field theories of the form Eq. [T4"|) were 
proposed at the quantum level and h appears explicitly; we are focused on the 
classical, h = 0, level. Note also that the scale of the coupling in Eq. ([TB")) is 
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set by £de- This is a macroscopic in length scale, and if we expand R about 
4:A DE G/c 2 , Eq. ([T3J reduces to Eq. (0 with a 

f = J)'(4)^2£, (15) 

which has a magnitude ~ 10 75 — indicative of an inherently classical theory — if 
</> has the mass of a proton. This value for £ is seventy-five orders of magnitude 
larger than the values of £ usually considered. It also signifies that a perturba- 
tive solution of Eq. (|13p would be of limited use at best, and the non-linearity 
of the coupling must be explicitly taken into account. 

2.4.2 The f(R) Theory 

Proposals for modifying the Hilbert action by considering functions, f(R), of 
the Ricci scalar have been made before (see [23] and [25J for reviews). These 
theories were first introduced to explain cosmic acceleration without the need 
for Dark Energy using a 1/R action [25], [27], and further extensions of this 
model have been made [25], [25] since then. They are now being studied in 
their own right, and various functional forms for f(R) are being considered. 
Indeed, connection to Modified Newtonian Dynamics has been made for log- 
arithmic f(R) terms [3D], [31], while with other choices of f(R) connection 
with quintessence has been made [32J - [37J as well. Importantly, issues with 
the introduction of a "fifth force" , and compatibility with terrestrial exper- 
iments have begun to be addressed through the Chameleon Effect (see [35] - 
[41] and an overview in [32] ) , which is used to hide the effects of field with a 
small mass that would otherwise be seen. 

It is also important to note that while f(R) theories change the action for 
gravity, in our approach we do not; we still take the action for gravity to be 
the Hilbert action with the addition of a cosmological constant. We instead 
change the response of matter to gravity by extending the equations of motion 
for test particles, and thus change the energy-momentum tensor for matter. 



3 The Energy-Momentum Tensor 

Beginning with [6j, there have been a number of attempts to show that the 
GEOM are a necessary consequence of the Einstein's field equations, Eq. (TTJ). 
Modern attempts at demonstrating such a linkage [7J , [5J focus on the energy- 
momentum tensor, and consider the motion of a test particle moving in a 
region of spacetime where T^ = outside of a "worldtube" that surround 
the test particle; inside this worldtube, Einstein tensor G^ v ^ 0. In fact, this 
tensor must satisfy the strong energy condition G^i'V < (for our signature 
for the metric) there, where < M and t' are two arbitrary, time-like vectors. As 
shown in [5J, the test particle then necessarily moves along a geodesic. While 
this proof do not explicitly include the cosmological constant term, replacing 
G>, y = G^ v + AdeQuv does not materially change the nature of the proof given 
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in [5]; since A E e > 0, Gpv satisfies the strong energy condition as long as G M „ 
does. 

If T^ v = then 9^[4 + %ttT/A de ] = SH[4], and is a constant. We then see 
explicitly from Eq. © that the extended GEOM reduces to the GEOM. As 
the T M „ = case is precisely the situation covered by [5] , the extended GEOM 
does not violate these theorems. Indeed, in the following we will explicitly 
construct the energy-momentum tensor for dust within the extended GEOM 
framework. 

Consider a collection of massive particles that can be treated as an inviscid 
fluid with density p, pressure p, and fluid velocity, v fi (x). We consider the 
spacetime to be spatially symmetric, so that the most general form that the 
energy-momentum tensor for this fluid is the usual 

T M „ = pv^Vu - (g M v ^J P- (16) 

We emphasize that this form for T^ v depends only on the spatial isotropy of 
the fluid, and thus holds for both the GEOM and the extended GEOM. 

Following |21j . energy and momentum conservation, VT^,, = 0, requires 
that 

= v v V u {p + p/c 2 )v,j, + (p + p/c 2 )V v v u v„ + (p + p/JWVvVp - V pP . (17) 

Since v 2 = constant even within the extended GEOM formulation, projecting 
the above along v^ gives once again the first law of thermodynamics 

d(Vpc 2 ) = -pdV, (18) 

where V is the volume of the fluid. This analysis holds for both the GEOM and 
the extended GEOM, and thus the first law of thermodynamics holds for both 
equations of motion. The standard analysis of the evolution of the universe 
under the extended GEOM therefore follows much in the same way as before. 
Next, projecting Eq. (fT?) along the subspace perpendicular to v^ gives the 
rclativistic version of Euler's equation 



= 



(p+5)/V^-( V -^)w (19) 



Once again, Eq. (|T9")) holds for both sets of equations of motion. 

Consider now the simplest case when the constituent test particles in the 
fluid do not interact with one another except under gravity. This corresponds 
to the case of "dust" . If test particles in this dust follow the GEOM, the 
solution to Eq. flU]) gives the usual T^ EOM - Dust = pv^Vu with p = 0. If, on 
the other hand, the test particles follow the extended GEOM, the situation 
changes. Using Eq. ((?]), Eq. (fT!?]) becomes 

o = (g^ - ^n) { ( P c 2 + P) V log n - V>} , (20) 

so that 

{pc 2 + P )V ti D\-mV ti p = Z DE A DE cv^ (21) 
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where (;de is a constant. By contracting the above with v^, it is straightforward 
to see that if £de 7^ 0, p will increase linearly with the proper time. This would 
be unphysical, and we conclude that £de must be zero. 

Taking the pressure as a function of only the density, Eq. ([2~Tj) reduces 
to a nonlinear, first order, ordinary differential equation. We will not solve 
this equation in general. Instead, we look at the nonrelativistic limit where 
pc 2 3> 3p. Then T sw pc 2 , and *H can be approximated as a function of p only. 
Equation (j2"i"T) can then be solved implicitly to give 

ff rig 

p(p) = -pc 2 + c 2 <K[4 + 8wp/A DE ] / . (22) 

J $H[4 + 8irs/A DE \ 

Given the density, the pressure is then determined once a form for 9\ is known. 
The energy- momentum tensor for dust, T E ' xt ~ I)ust , under the extended GEOM 
can then be constructed using p, and the resultant p from Eq. I|22|l . Note 
also from Eq. ([22]) that p « 47r^'[4]p 2 c 2 /9t[4]A DE when p ->■ 0, while p « 
-yl DS c 2 (1 - 1/«H[4]) when p > /b B /27r. Thus, T^- dust « /ro^, and the 
solution Eq. (j2"2")) is consistent with the approximation T f» pc 2 . 

The physical reason for the presence of this pressure term in j^xt-Dust can 
be seen from Eqs. (J7J) and IjHJ). If the second equation is solved for the energy 
of a particle, it is straightforward to see that under the extended GEOM a 
collection of test particles behave as though they were in an external potential 
set by i D(c 2 R/AdeG)- As such, the particles no longer follow geodesies as 
they do in the GEOM even though they only interact with each other through 
gravity. The presence of the nonzero pressure term in yE^t-Dust j g a renec ti on 
of the presence of this potential. 



4 A Form for 2) (x) and Experimental Bounds on a. a 

Our analysis up to now is valid for all S(x). Requiring that our extension of 
the GEOM does not produce effects that would have already been observed 
in experiments will fix a specific form for 2)(x). 

Since our extension of the GEOM does not change the equations of mo- 
tion for massless test particles, we expect Eq. (JTJ) to reduce to the GEOM 
in the ultrarelativistic limit. It is thus only in the nonrelativistic limit where 
the effects of the deviations from the GEOM due to the additional terms in 
Eq. ||ZJ) can be seen. We therefore focus on the impact of the extension in the 
nonrelativistic, weak gravity limit, and begin by expressing Eq. (J7J in these 
limits. 



4.1 Constructing J)(x) 

We first perturb off the Minkowski metric r/^ in the weak gravity limit by 
taking g^ v = r/^ + /i M „, where the only nonzero component of h^ v is ft-oo = 
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2<£/c 2 , and ^ is the gravitational potential. Equation ([T]) then gives 

V 2 $ + 2^^-$ = 4TrpG-A DE G, (23) 

in the presence of a cosmological constant. Next, the temporal coordinate, x , 
for the extended GEOM in this limit will, as usual, be approximated by ct to 
lowest order in |v|/c. The spatial coordinates, x, on the other hand, reduce to 



cPx _ _ / Attc 



V<P 



dt 2 \A 



IDE 



®'(4 + 8irp/A DE ) 



1 + D (4 + 8np/A DE ) 



Vp. (24) 



Here, we have assumed that the spacetime is spatially symmetric, and that 
the particle moves through an ambient, nonrelativistic fluid with density p. 

For the additional terms in Eq. (f2"4"|) from the extension not to contribute 
significantly to Newtonian gravity under current experimental conditions, S'(4-t 
8ttp/Ade) -> when p » Ad E /2it. We can now either choose £>'(ir) > or 
'X)'(x) < 0, as 1)'(x) — > from either above or below. This is an ambiguity in 
the construction of T)(x) that cannot be resolved with our current analysis; 
for now, we take 1)'(x) < 0. The simplest form for T)'(x) with the correct 
asymptotic properties is 

V'{x) = -- \^, (25) 

where a a is a power-law exponent that determines how fast 1)'{x) — > for 
x — > oo, and x is t ne normalization constant 



f f 00 ds 



X Jo 1 + s 1+aA ' 



(26) 



To ensure nonzero effective masses, 2) (x) must be positive, and we integrate 
Eq. ([25]) to get 

with the condition that a a > to ensure that the the integral is defined. With 
this choice of integration constants, 3D(0) = 1 and S(s) ->0asa;^oo. While 
the precise form of 25 (x) is calculable, we will not need it. Instead, because 
8irp/A DE > 0, 

T>(4 + 8np/A DE ) = x J2 JTZ YZ 4 + T^ ' ( 28 ) 

^ ?i(I + a A ) + ola \ A DE 

while 

I = 1 + 2 V C -1 )" (29) 

X ^ [I + (n+l)(l+ aA )][n(l + ^) + ^]- l } 

Notice that in the a a — >■ oo limit, 33 (a;) — > 0, £Ext — > £o, and the GEOM is 
recovered. 
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Bounds on a a will be found below. For now, we note that for a a > 1, 

X ~ 1 and D (4 + 8irp/A DE ) w 0. Thus, 



1+omI j 



£~-™ + ra i+ «+& < - <3o) 



4.2 Constraints From Terrestrial Experiments 



We now consider the constraints placed on 1)(x) due to terrestrial experiments 

1.82 
-0.84 



From WMAP, A DE = 7.21°§ 2 84 x 10 -30 g/cm 3 , which for hydrogen atoms 



corresponds to a number density of ~ 4 atoms/m 3 . It is clear that the den- 
sity of both solids and liquids far exceed A EE , and in such media Eq. (f5U|) 
reduces to what one expects for Newtonian gravity. Only very rare gases, in 
correspondingly hard vacuums, can have a density that is small enough for 
the additional terms in Eq. (J24I) to be relevant. The hardest vacuum currently 
attainable experimentally is ~ 1CP 13 torr [42]. For a gas of He4 atoms at 4 
°K, this corresponds to a density of pumit ~ 10 -18 g/cm 3 , which is 11 orders 
of magnitude smaller than Aoe- Nevertheless, because the scale of the ac- 
celeration from the additional terms in Eq. (150)1 — which goes as c 2 /L for the 
relevant length scale L in the experiment — is so large, effects at these densities 
will nevertheless be relevant. Indeed, we will use pnmit to determine a lower 
bound to ola- 

Consider a simple experiment that looks for signatures of the extension of 
the GEOM Eq. (|30l) by looking for anomalous accelerations (through pressure 
fluctuations) in a gas of He 4 atoms at 4 °K with a density /Oiimit- Inside this gas 
we consider a sound wave with amplitude epnmit propagating with a wavenum- 
ber k. Suppose that the smallest measurable acceleration for a test particle in 
this gas is abound- For the additional terms in Eq. (J5U)) to be undetectable, 



A ( . Al ' l: 



a a 



abound > ~PT o ek - ( 31 ) 

2 \ 87rpii m it / 

This gives a lower bound on a a as 

_ log [2a bo und/c 2 xefc] , . 

a /lbound — 1 fl Jo T- \ 6Z > 

log [A DE /8-iTpii m it\ 

For e = 0.1, k = 1 cm , and abound = 1 cm/s 2 , a/ibound ranges from 1.28 
for A DE = 10 -32 g/cm 3 to 1.58 for A DE = 1CT 29 g/cm 3 . Notice that because 
a/ibound depends only logarithmically on e, fc, and abound, the lower bound for 
a a is relatively insensitive to the specific values taken for these parameters. 

Adelberger has recently done state-of-the art, Eovtos-type experiments 
[43] . [44], [45] to test the Newtonian, 1/r 2 law for gravity. While the pressures 
under which these experiments were performed were not explicitly stated, as 
far as we know these experiments were not done at pressures lower than 10 -13 
torr; we thus would not expect effects from extended GEOM to be apparent 
in these experiments either. 
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4.3 Constraints from Solar System Observations 

Most of the classical tests of general relativity are based on the motion of 
bodies in the solar system. We now consider the constraints placed on 2)(x) 
due to these tests. 

We begin by noting that in idealized situations such as the analysis of 
the advancement of perihelion of Mercury, the energy-momentum tensor is 
taken to be zero outside of a massive body such as the Sun; the terms on 
the right hand side of Eq. (0 will clearly not affect these analyses. While 
this argument would seem to hold for this and the other experimental tests of 
general relativity, it is much too simplistic. In practice, the T^ in each of these 
tests does not, in fact, vanish; a background density is always present. Except 
for experiments involving electromagnetic waves, what is instead needed is a 
comparison of the background density with Aee- It is only when this density 
is much greater than Aee/^tt that the additional terms in Eq. (0 may be 
negligible. We present this analysis below. 

Of the classical, solar system tests of general relativity, only in observa- 
tions or experiments involving motion of massive test particles can effects of 
the extension be seen. The extended GEOM does not affect the motion of 
electromagnetic waves, and such tests of general relativity as the gravitational 
redshift and the deflection of light by the Sun will not be affected by the ex- 
tension. Considering first the advancement of the perihelion of Mercury, we 
note that the density of matter, p rb, (due primarily to the solar wind) in the 
region of Mercury's orbit is ~ 10 -23 g/cm 3 . This density is only five orders 
of magnitude smaller to pumit- There is, then, a possibility that the extended 
GEOM could change the amount in which the perihelion of Mercury advances 
due to general relativistic effects. 

To show that our extension nevertheless does not change the advance of 
perihelion, we follow extensively the derivation of the advancement of perihe- 
lion given in [3T]. In particular, we use a Hamilton- Jacobi-based analysis to 
calculate the advancement and replace Eq. © by 



t M t 



1 + 2) 4 



8irp, 



'orb 



A 



DE 



GM t n (GM 



\ c 2 



dSi 



M 



cdt 



GMn 



o 



'dS 



M 



Or 



ldS } 



M 



r d(f> 



(33) 



Here, Tom is the mass of Mercury, Mq is the mass of the Sun, and Syi is 
the action for the motion of Mercury about the Sun. We have used the usual 
Schwarzchild solution for the region outside of the Sun, and kept terms up 
to second order in GM/c 2 r. While in principal we would expect the extended 
GEOM to change the solution of Einstein's equations, we are working in both 
the weak gravity and nonrelativistic limits, where, as we have shown above, 
Port>c 2 ; the Schwarzchild solution is unchanged in these limits. We 



T 



Ext-Dust 
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also note that because p rb ~ 10 -23 g/cm 3 , the additional matter in Mercury's 
orbit will not affect the form of the Schwarzchild solution, and can be neglected 
when determining the metric outside of the Sun. 
We next look for solutions of the following form: 



Sm = -m M c 2 Et + m M Lip + m M c &(r), 



(34) 



where E is the energy per unit rest energy, L is the angular momentum per 
unit rest mass, ip is the azimuthal angle, and r is the radial position of Mercury 
measured from the Sun. Solution of Eq. (|3"3")l then gives 



6 = ± 



CM r 
(l + £-£ 2 ) + 2-^- l-2(l-E 2 ) 



L 2 



(GM\ 



1/2 



6 



dr. 



(35) 



after keeping terms up to order (GMq/c 2 t) 2 . In obtaining Eq. (1351) . we have 
used the fact that 1 — E 2 ~ GA/q/c 2 o, where a is the semi-major axis 
of Mercury's orbit; observations set a — 57.91 x 10 6 km [46], giving 1 — 
E 2 ~ 1CT 8 . We next used the relation from Eq. (gSJ), 25[4 + 8iTp olh /A DE } « 
(x/oca) (A D E/&Trporb) aA - Thus, £> - lO" 10 - 1CT 12 for a A = 1.28 - 1.58, so 
that GM Q /c 2 a < T)< (GM Q /c 2 a) 2 . 

The shape of the orbit of Mercury is determined by minimizing the action, 
dS/dL = 0, in Eq. ((34)) . The resultant integral is straightforward to calculate, 
and we obtain 



<p = 



1 



( GM Q 
V cL 



cos 



(l-e 2 )a 1 
e 



er 



(36) 



where analytically 

_ GM Q 1 - 2(1 - E 2 ) 



cL 



1 + 25 

\ 2 



GMr. 



O 



E 2 
1 + 25 -E 2 



i + 4(i-£; 2 )-6 



-J 



GM, 



o 



V ci 



(37) 



and e is the eccentricity of the orbit. The advancement in the perihelion of 
Mercury therefore still has the form 



Sip = 



GttGM, 



O 



c 2 a(l — e 2 ) 



(38) 



Note, however, that the product a(l — e 2 ) is independent of 25. Thus, the 
amount that the perihelion of Mercury 's orbit advances due to general rela- 
tivistic effects is not changed by the extended GEOM. Physically, this is be- 
cause the extended GEOM only modifies the rest mass in Eq. (0, and does 
not modify terms explicitly dependent on the velocity of Mercury. 
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What the extended GEOM does affect is the analyticai expressions for 
e and a. However, in calculating the numerical value of Sip, both e and a 
are taken as measured quantities obtained from observations; they are not 
calculated from first principles. The fact that there is now a slight different 
relationship between e and a, and the total energy, mMC 2 E, and the angular 
momentum, itlmL, of Mercury (by less than 0.01 % for a a = 1.58) would 
require an independent method of determining E and L to check. Such an 
independent determination is not currently available. 

There are three other potential solar system tests of general relativity that 
could, in principal, be used to constrain the extended GEOM. The first is 
looking at three-body effects in the lunar orbit due to the Sun, Earth and the 
Moon. It is questionable whether such effects can be measured [31], however, 
and they therefore cannot serve as a check on the extended GEOM. The second 
is frame dragging, where the precession of a spinning object is caused azimuthal 
changes in the local metric of spacetime by the rotating Earth; this effect 
is being measured by Gravity Probe B [47]. The third also measures frame 
dragging, but uses two satellites — LAGEOS and LAGEOS 2 — instead of a 
spinning body, to detect the Lense-Thirring effect; frame dragging was recently 
seen [35] in this experiment. 

Frame dragging, whether it is by a spinning object or by two orbiting satel- 
lites, is inherently a velocity-dependent effect, however, that couples either to 
the spin of the object or to the orbital angular momentum of the satellites 
with the spacetime metric. Our extension of the GEOM changes the rest mass 
of the test particle, and in a spherically symmetric geometry in the nonrel- 
ativistic limit, the additional effects due to the extension is radial. We thus 
would not expect the effects from our extension can be seen either from the 
precession of a spin, or through the Lense-Thirring effect. 



5 Concluding Remarks 

We have shown that because of the existence of a universal length scale, £de, 
it is now possible to construct an extension of the GEOM. This extension 
preserves the equivalence principal, does not change the motion of massless 
test particles, and does not produce effects that would be detectable in either 
terrestrial experiments, or through observing the motion of bodies in the solar 
system. Our extension of the GEOM is thus a physically viable alternative to 
the GEOM. 

The question remains as to whether these equations of motion have any 
physical relevance. In short, is anything gained by using this extension? Be- 
cause £de = 14010|2o Mpc, we would expect that any effects from the ex- 
tended GEOM will become apparent at much longer length scales than those 
considered here. Indeed, given the size of Ide the only reason why we would 
expect the extended GEOM to be relevant at all is because 2) is a nonlinear 
function of the energy-momentum tensor of ambient matter. This question of 
relevance will be addressed in a future paper [15] where the extended GEOM 
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is applied to the motion of stars in the rotation curves of galaxies, and to 
the density of matter at cosmological length scales. These are the scales at 
which we expect the effects from the extension to come into play, and where 
its relevance can be assessed. 

Finally, as noted in [21] , Eq. (fl9)l can be solved in general to give an equa- 
tion of motion for particles, v v V u Vfj, — V IM p/(p + p/c 2 ), and we see that the 
presence of any pressure term in the energy-momentum tensor results in de- 
viations from geodesic motion. Given that Ajje can also be used to construct 
a pressure, it is natural to ask whether the effects of the extended GEOM 
can be obtained through the introduction of an ad hoc pressure term in the 
energy-momentum tensor. Such an ad hoc term can only be introduced to the 
energy-momentum tensor for matter, however; for the reasons given in the in- 
troduction, the equations of motion for massless particles cannot be changed. 
In addition, because the behavior of any massive particle approaches that of 
a massless one in the ultrarelativistic limit, this ad hoc pressure term must 
be constructed so that irrespective of frame this term contributes a negligible 
amount to the energy of the particle in this limit. Moreover, even if such a con- 
struction can be accomplished at the ultrarelativistic limit, hurdles remain at 
the nonrelativistic limit. While it is possible to construct in the nonrelativistic 
limit an appropriate ad hoc pressure using Ade and 2), when the resultant 
equations of motion are applied to the same systems as the extended GEOM 
in \T5\ , effects are predicted at the cosmological scale that either do not agree 
with experiment or are not physically reasonable. This occurs even though the 
pressure is chosen so that at galactic scales predicted effects will be in broad 
agreement with observations. For these reasons, it is doubtful that introducing 
ad hoc pressure term in place of the extended GEOM will be succussful. 
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